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Abstract 

For two massless particles i and j. the collinear limit is a special kinematic configuration 
in which the particles propagate with parallel four-momentum vectors, with the total mo¬ 
mentum P distributed as pi — xP and pj — (1 —x)P, so that Sjj = ( pi + Pj ) 2 — P 2 — 0. 
In Yang-Mills theory, if i and j are among N gauge bosons participating in a scattering 
process, it is well known that the partial amplitudes associated to the (single trace) group 
factors with adjacent i and j are singular in the collinear limit and factorize at the lead¬ 
ing order into (IV—l)-particle amplitudes times the universal, x-dependent Altarelli-Parisi 
factors. We give a precise definition of the collinear limit and show that at the tree level, 
the subleading, non-singular terms are related to the amplitudes with a single graviton 
inserted instead of two collinear gauge bosons. To that end, we argue that in one-graviton 
Einstein-Yang-Mills amplitudes, the graviton with momentum P can be replaced by a pair 
of collinear gauge bosons carrying arbitrary momentum fractions xP and (1 —x)P. 


In high energy particle physics, collinear kinematics are very common. Viewed in the 
laboratory frame, all quarks and gluons (partons) propagating inside protons accelerated 
at the Large Hadron Collider (LHC) move in the beam direction, with a very little of 
transverse momentum. Since the early days of Quantum Chromodynamics (QCD), such 
collinear parton configurations have been in the focus of perturbative computations. In 
the so-called leading logarithmic approximation, the violation of Bjorken scaling [I] in 
proton structure functions can be understood as an effect of 1 —* 2 parton decays which 
are necessarily collinear. They are described by Altarelli-Parisi probabilities [2] and involve 
the running gauge coupling constant that brings the fundamental QCD mass scale. When 
protons collide at high energies, many quarks and gluons are often produced in a single 
two-parton collision. Multi-parton amplitudes favor collinear final state configurations due 
to the singular behaviour that will be discussed below. At the LHC, such partons fragment 
into hadronic jets. 

In general, for two massless particles i and j, the collinear limit is defined as a special 
kinematic configuration in which the particles propagate with parallel four-momentum 
vectors, with the total momentum P distributed as pi — xP and pj — (1 — x)P , so that 
Sij = (pi + Pj ) 2 — P 2 = 0. In QCD, as in any Yang-Mills theory, if i and j are among N 
gluons participating in a scattering process, it is well known that the partial amplitudes 
[3] associated to the (single trace) color factors with adjacent i and j are singular in 
the collinear limit and factorize at the leading order into (N— l)-gluon amplitudes times 
the universal, .x-dependeni Altarelli-Parisi factors (three-gluon MHV amplitudes). The 
singularity is a simple pole at .sy ? — 0 due to an intermediate gluon propagating on zero 
mass shell. In this paper, we go beyond the leading pole approximation and discuss non- 
factorizable, finite contributions 1 . 

In order to give a precise definition of the leading and subleading parts, we need to 
specify how the collinear limit is reached from a generic kinematic configuration. Let us 
specify to generic light-like momenta pi — Pn-i, Pj — Pn and introduce two light-like 
vectors P and r such that the momentum spinors decompose as 

A n -i = A p cos 9 — eA r sin 9 , A n-i = Ap cos 9 — e\ r sin 9 , 

Ajv = Ap sin 9 + e\ r cos 9 , A at — Xp sin 9 + eX r cos 9 , 


hence 


Pn-i = c 2 P — esc(XpX r + X r Xp) + e 2 s 2 r , 
Pn = s 2 P + esc(ApA r + A r Ap) + e 2 c 2 r , 


( 2 ) 


1 The “leading logarithms” come from integrating such poles. In the language of perturbative 
QCD, the subleading terms discussed here belong to so-called “higher twist” contributions. 
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where 

c = cos 6 = y/x , s = sin 6 — yjl — x . 

We also have 

(N-l N) = e ( Pr) , [N-l N] = e [Pr] . 
The total momentum is: 


(3) 

(4) 


Pn-i +PN = P + e 2 r , (pn-i +Pn ) 2 = Sjv-i,jv = 2Pre 2 . (5) 


The collinear configuration will be reached in the e —> 0 limit and the tree amplitudes 
discussed below will be expanded in powers of e. Partial amplitudes with adjacent N—l 
and N contain single (factorization) poles. Thus their leading terms are of order 0(e~ 1 ) 
and the subleading ones are of order O(e 0 ). Collinear expansions of partial amplitudes 
with non-adjacent N—l and N start at the subleading O(e 0 ) order. 

The leading collinear behaviour of amplitudes with adjacent N—l, N is well known [3] 
and depends on respective helicities. For identical helicities, 

n(l,..., JV-1 + , JV+) = — 1 A(l,.. . , P+) + e°Sub ++ + . .. , 

(N—1N) sc (6) 

A(l,...,N-l-,N-)= |JV _ ljv]s / (l.-".-r-) + e°Sub + ... , 

where we used superscripts to denote helicity states. Here, Sub denote subleading contri¬ 
butions which are the focus of this work. The remaining terms vanish in the e —>• 0 limit. 
For opposite helicities: 


A(l, ■ ■ ■, N—l + , N~) — 


(N-lN)c 


A(1,...,P~) + 


[iV-liV]s A ^ 1, '" ,jP+ ' 1 (7) 


+ e u Sub 


+- 


+ ... 


The starting point for our discussion of subleading terms is the recent observation [4] 
that the tree-level Einstein-Yang-Mills amplitudes describing decays of a single graviton 
or a dilaton into a number of gauge bosons, can be written as linear combinations of pure 
gauge amplitudes in which the graviton (or dilaton) is replaced by a pair of gauge bosons. 
Their ±1 helicities add up to ±2 for the graviton or to 0 for the dilaton. Each of them 
carries exactly one half of the original graviton or dilaton momentum, which is a special 
case of a collinear configuration with s = c = yj 1/2. From now on we will focus on 
graviton amplitudes. The crucial point is that the relations derived in [4] can be extended 
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to arbitrary collinear configurations, in the following way 


-Aeym(I) 2,..., N— 2; P ±2 ) = ( 8 ) 

2 rrfl -1 I I 

51 l Z =2 *=2 j=i 

N -3 N—2 i 

+ 22 ^2 ( ^ym(1, • • • > l, N— l^, l + 1,..., i, N ± , i + 1,..., N—2) >, 

where n and g are the gravitational and gauge coupling constants, respectively 2 . On the 
left hand side, we have a mixed gauge-gravitational amplitude involving a single graviton 
of momentum P, helicity ±2 as indicated by the superscript, and N—2 gluons. This 
amplitude is associated to a single trace color factor with the respective gluon ordering. 
On the right hand side, we have a linear combination of pure gauge, partial amplitudes 
weighted by the kinematic invariants s^jv-i = ZpjPN-i- Here, the graviton is replaced by 
two gluons in the collinear configuration: 

Pat-i = c 2 P = xP , = s 2 P = (1 — x)P , (9) 


i.e. the leading O(e 0 ) order of Eqs. (1) and (2). Note that on the right hand side, N —1 
and N are never adjacent, therefore the Einstein-Yang-Mills amplitude emerges from the 
collinear limit of Yang-Mills amplitudes at the subleading O(e 0 ) order. In order to further 
discuss Eq. ( 8 ), it is useful to write it explicitly for — 5, 6 , 7: 


H(1,2,3;P ±2 ) - ^ s 24 71(1,5*2,4^3) , 

r 

2 

H(1,2,3,4;P ±2 ) = ^-|s 25 4(1,6*, 2, 5*, 3,4) + S 45 4(1, 2, 3, 5*, 4, 6 *)} , 

2 

4(1, 2, 3,4, 5; P ±2 ) - {s 26 A{ 1, 7*, 2 , 6 *, 3,4, 5) + s 36 A(l, 2 , 7*, 3, 6 *, 4, 5) 

+ (s 3 e + « 2 e) A( 1 , 7*, 2 , 3, 6 *, 4, 5) + s 56 A( 1, 2 , 3, 4, 6 *, 5, 7*) | 


( 10 ) 

( 11 ) 

( 12 ) 


The fact that the relations written in Ref. [4] can be extended from s = c = \J 1 / 2 , 
i.e. from x — 1 / 2 , to an arbitrary collinear configuration by inserting a simple factor 
of s 2 = 1 — x is highly non-trivial. It is easiest to check for the helicity configurations 
described by MHV amplitudes [5] on the r.li.s. of Eq. ( 8 ), i.e. when the collinear pair is 

2 is the smallest integer greater than or equal to . Since the graviton is identified by its 

momentum P, we can skip in the following the EYM and YM labelings of the amplitudes. 
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among TV — 2 gluons with identical lielicities and there are only two gluons with opposite 
helicities 3 . Then the r.h.s. of Eq. (8) is a homogenous function of spinor (and momentum) 
variables and it is easy to see that, in this case, arbitrary value of x can be reached 
from a; = 1/2 by a simple rescaling of the amplitudes, with the net effect of an overall 
1 — x factor. For other helicity configurations, the amplitudes are not homogenous in the 
momenta of collinear gluons. Already at the NMHV level, individual amplitudes contain 
poles in three-gluon channels [ijN— 1] and [ij TV] (with i. j ^ TV—1,TV), characterized by 
the kinematic invariants 


tijN—l — (Pi T Pj T PN—l) C tijp T S Sij T ... 

/ \2 2 2 \ ^/ 
tijN = (Pi + Pj + PN) = S tijp + C Sij + . . . 

Such poles must cancel on the r.h.s. of Eq. ( 8 ) for the EYM amplitude to be free of 
unphysical singularities. In the appendix, we show that it is indeed the case for TV = 6 , 
and obtain an explicit expression for A(l + , 2 + , 3“, 4 + ; P~ 2 ) in agreement with Eq. (11). 
For IV = 7, a similar check is still possible but it involves very tedious computations. 
Starting from TV — 8 , NNMHV amplitudes can appear on the r.h.s. of Eq. ( 8 ), therefore a 
complete proof would have to rely on more general representation of tree amplitudes or on 
recursion relations. Actually, the most straightforward way is to consider these amplitudes 
as a zero-slope limit of superstring disk amplitudes involving open and closed strings. 
Then, Eqs. ( 8 ) can be proven for arbitary helicity configurations [ 6 ]. In this work however, 
we focus on field-theoretical amplitudes. 

At the tree level, there are (TV — 3)! independent TV-gluon amplitudes [7]. For a 
given TV, we can express the Yang-Mills amplitudes appearing in Eq. ( 8 ) in terms of the 
basis A( 1, cr(2, 3,.. .IV—2), TV—1, TV), where a denotes the set of (TV — 3)! permutations of 
2, 3,..., TV—2. Let us start from TV = 5, as in Eq. (10), where we use: 

A(l,5,2,4,3) = A(l,2,3,4,5)+ g21 + g23 A(l,3,2,4,5) . (14) 

-§25 s 25 

As a result, we conclude that the following relation is valid up to the O(e 0 ) order: 

s 3P A( 1, 2 , 3, 4±, 5±) - s 2P A(1, 3, 2 , 4±, 5±) = ^-A( 1, 2 , 3; P ± 2 ) . (15) 

KX 

By using Eq. ( 6 ) and BCJ relations for four-gluon amplitudes [7], it is easy to see that the 
leading collinear singularities 0{e~ l ) drop out, therefore Eq. (15) connects the subleading 

3 The other case, when the collinear pair carry helicities opposite to all other TV—2 gluons, 
does not contribute because the corresponding amplitudes vanish in the collinear limit as e 4 . 
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terms with the mixed gauge-gravitational amplitude. For N — 5, we obtain one relation 
between the subleading parts of two independent amplitudes. For N == 6 , a similar equation 
reads: 

S4pA(1, 2, 3, 4, 5 ± ,6 ± ) - s 3 p[dl(l,2,4,3,5 ± ,6 ± ) + A(l,4,2,3,5 ± ,6 ± )] 

a 2 (16) 
+ S 2P A(1, 4, 3, 2, 5 ± , 6 ± ) = —A( 1, 2, 3,4; P ±2 ) . 

KX 

In this case, however, we have two additional mixed amplitudes, say A{ 1, 3, 2,4; P ±2 ) and 
A(l, 2, 4, 3; P ±2 ), that can be used in similar relations, obtained by interchanging 2 -H- 3 
and 3f>4, respectively. As a result, we obtain three relations for the subleading parts of 
six independent gauge amplitudes 4 . For N — 7, 


s 5P A( 1,2,3,4,5,6 ± ,7 ± ) 

- s 4P [A(1, 2 , 3, 5,4, 6 ± , 7 ± ) + A( 1 , 2 , 5, 3,4, 6 ± , 7 ± ) + A(l, 5, 2 , 3, 4, 6 ± , 7 ± )] 

+ 8 3P [A(1, 5, 4, 2 , 3, 6 ± , 7 ± ) + A( 1 , 5, 2 ,4, 3, 6 ± , 7 ± ) + A(l, 2 , 5, 4, 3, 6 ± , 7 ± )] (17) 

a 2 

- s 2 pA(l,5,4,3,2,6 ± ,7 ± ) = —A(l, 2, 3, 4, 5; P ±2 ) . 

KX 


In this case, there are 24 independent Yang-Mills amplitudes with the subleading collinear 
behaviour constrained by twelve Einstein-Yang-Mills amplitudes. For arbitrary N a similar 
formula reads 


E (-l) mp s p(N . 2 )P A(l, p( 2, ..., N - 2), N - 1, N) 

P&Vn 



KX 


A(l,..., N — 2; P) , (18) 


where "Pat is a subset of permutations acting on 2,..., N — 2 and m p G {0,1} as specified 
in [ 6 ]. Now there are (N — 3)!/2 independent constraints. 

We see that the subleading collinear behaviour of pure gauge amplitudes is determined 
in part by the amplitudes with the graviton inserted instead of the collinear pair. Twice as 
many constraints are necessary, however, in order to fully determine the subleading terms 
for all amplitudes. In another physically interesting case of soft (x —> 0) divergences, the 
subleading behaviour has been recently discussed in Einstein’s gravity [ 8 ] and in Yang- 
Mills theory [9]. We hope that similar considerations will allow complete determination of 
the subleading behaviour in the collinear case. 

The fact that the graviton can be replaced by two gluons in arbitrary collinear config¬ 
urations in the single-graviton amplitudes of Eq. ( 8 ) raises an interesting question whether 

4 Three other mixed amplitudes are related by parity reflections, therefore they do not provide 
additional constraints. 
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pure Einstein, multi-graviton amplitudes share this property. The recent linearization [10] 
of KLT relations [11] suggests that this may be the case. It would be another indication 
for the existence of some underlying gauge structure in quantum gravity. 


Appendix 

We will show that Eq. (11) holds for A(l + , 2 + , 3“, 4 + ; P~ 2 ). To that end, we take the 
collinear limits, c.f. Eq. (1), of the six-gluon NMHV amplitudes written in Ref. [3]: 


A(l+,6-,2+,5-,3-,4+) 


A(l+,2+,3-,5-,4+,6-) 


(3P) 4 J [14] 2 [23] 2 | [12] 2 [34] 2 

S 1 PS 2 PS 3 P ) C 2 (S 2 S 2 3 + C 2 Si 4 )si4 S 2 (s 2 S 34 + C 2 S 42 )s 34 

+ [12] [23] [34] [41] ) | 
c 2 s 2 s 44 s 34 J 

(19) 

(3P) 4 / [14] 2 [23] 2 [12] 2 [34] 2 

SipS3pS 4P \ S 2 (S 2 S 2 3 + C 2 S 44 )S 23 C 2 (S 2 S 34 + C 2 S 12 )si 2 


[12] [23] [34] [41] 

C 2 S 2 S 42 S 23 

( 20 ) 

where we omitted terms that vanish in the e —> 0 limit; we also set g 2 = n — 1. After 
substituting into Eq. (11) and using momentum conservation, we obtain 



A(l + ,2 + ,3“,4+;P- 2 ) 


(3P) 4 

(12)(23)(34)(41) ’ 


( 21 ) 


in agreement with Ref. [4], For other NMHV helicity configurations, Eq. (11) follows in 
exactly the same way. 
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